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\ Abstract. We derive eigenvalue asymptotics for Sturm-Liouville operators with sin- 

fSJ ■ gular complex-valued potentials from the space W.f _1 (0, 1), a e [0, 1], and Dirichlet 

or Neumann-Dirichlet boundary conditions. We also give application of the obtained 
results to the inverse spectral problem of recovering the potential from these two 
\ spectra. 
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In this paper we shall study eigenvalue asymptotics for Sturm-Liouville operators 
on the interval [0, 1] with distributional potentials. Namely, we assume that q is a 
complex- valued distribution from the Sobolev space W 2 a_1 (0, 1), a G [0, 1], and consider 
an operator T that (formally) corresponds to the differential expression 

(1-1) Kf)-=-f" + qf 

and, say, Dirichlet boundary conditions. Explicitly, the operator T is defined by the 
regularisation method that was suggested in [4] for the particular potential q(x) = 1/x 
and was developed by Savchuk and Shkalikov in [30] for the class of distributional 
potentials from W 2 ~ (0, 1). (Incidentally, in this situation the form-sum [1, 10] and the 
generalized sum [20] methods yield the same operator). We observe that the considered 
class of singular potentials include Dirac 5- type and Coulomb 1/x- type interactions that 
are widely used in quantum mechanics and mathematical physics; see also [18] for other 
physical models leading to potentials from negative Sobolev spaces. 
The regularisation method consists in rewriting (|l.lj) as 

>< i (1.2) 1(f) = W) := -(/' - af)' - af, 

where a is any distributional primitive of q. We fix one such primitive a G L 2 (0, 1) 
in what follows and call the expression /' — af =: /M the quasi-derivative of the 
function /. The natural L 2 -domain of l a is 

dom/,, = {fe W}(0, l) | g Wt(0, l), Uf) e £2(0, l)}, 

and we observe that for / G doml^ the derivative /' = af + belongs to L 2 (0, 1) 
(but need not be continuous), so that dom/ CT C W^O, 1). 

In the present paper, we shall only focus on Sturm-Liouville operators T D = T a j) 
and Tn = T CT) n that are generated by l a and the Dirichlet and the Neumann-Dirichlet 
boundary conditions respectively, although other boundary conditions can also be 
treated in a similar manner (see, e.g., [18] for periodic and [31] for general regular 
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boundary conditions). In other words, Td and Tn are the restrictions of l a onto the 
domains 

domT D = {fe doml a | /(0) = /(l) = 0}, 
(L3) domT N = {/ e dom/ CT | /N(0) = /(l) = 0}. 

It is known [30] that the operators Td and Tn are closed, densely defined and have 
discrete spectra tending to +oo. We denote by X n (resp., the eigenvalues of Td 
(resp., Tn) counted with multiplicities and arranged by increasing of the real — and 
then, if equal, imaginary — parts of A„ (resp., /i n ). For definiteness, we shall always 
take A n and fi n from the set 

(1.4) VI := {z G C | -7r/2 < arg^ < vr/2} U {z = 0}. 

If a = 0, i.e., if a G L 2 , then the numbers A n and //„ obey the asymptotics [15,29-31] 

(1.5) A n = 7m + A n , /i n = 7r(n - |) + Ji n , 

where (A„) rag N and (jj, n ) n m are some ^-sequences. It is reasonable to expect that 
if a becomes smoother, then the remainders A n and Jl n decay faster; for instance, if 
a = 1, i.e., if q G £2(0, 1), then the classical result (see, e.g., [26, Theorem 3.4.1] or [28, 
Theorem 2.4]) states that A n ,/I n = O^n^ 1 ). Thus the problem arises to characterize 
the decay of A n , Ji n depending on a G [0, 1]. 

Our interest in the above problem has stemmed from the inverse spectral theory for 
Sturm-Liouville operators with singular potentials. Namely, we proved in [16] that, 
as soon as the numbers A^ and fj/ n £1X6 real, strictly increase with n, interlace, and 
obey (jl.5|) with £ 2 -sequences (A n ) and (Jin), then there exists a unique real- valued 
o G £2(0, 1) such that {A^} and {n n } are spectra of the Sturm-Liouville operators 
To-d and T^n respectively with distributional potential q = o' G W / 2 " 1 (0, 1). It is 
reasonable to believe that if A ra and Jt n decay faster, then a will be smoother. For 
example, the classical result of Marchenko [26, Theorem 3.4.1] claims that if, under the 
above assumptions, we have in addition 

n r\ \ A . ^ ~ A Ji' n 

(1-6) A n = 1 , /in = 1 

n n n n 

with real A and £ 2 -sequences (X' n ) and (p' n ), then the corresponding potential q is 
in L 2 (0, 1) (thus a = 1) and A = j- q, cf. also [24, Theorem 3.3.1]. It would be 
desirable to "interpolate" between a = 1 and a = and solve the inverse spectral prob- 
lem for all intermediate a G (0, 1). The essential step towards such project is to derive 
eigenvalue asymptotics for Sturm-Liouville operators with potentials in W / 2 o_1 (0, 1)— 
i.e., to treat the direct spectral problem. And indeed, based on the results obtained 
here, we completely solve the inverse spectral problem for Sturm-Liouville operators 
with potentials in the scale V4 / 2 Q_1 (0, 1), a G [0, 1], in our paper [17]. 

Another motivation for this work is the recent papers [18] and [31], where similar 
questions are considered. In particular, Kappeler and Mohr in [18] found eigenvalue 
asymptotics for the Dirichlet and periodic Sturm-Liouville operators with complex- 
valued potentials that are periodic distributions from the space W^^iO, 1), a G (0, 1]. 
The Dirichlet eigenvalues X n were proved there to obey the asymptotics 

(i.7) x i = „ v+m -!tMptM + „ m 
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where q{n) is the n-th Fourier coefficient of q and the sequence (u n ) belongs to l 2 ^ 16 
with e > arbitrary (the weighted t s v spaces are defined at the end of Introduction); 
see [18] for more precise formulations. The authors performed the Fourier transform 
to work in the weighted £2 spaces rather than in the Sobolev spaces and then derive 
the estimates for the resolvent that yielded the detailed localisation of the spectrum. 

Savchuk and Shkalikov [31] considered Sturm-Liouville operators with complex- 
valued potentials q that are distributional derivatives of functions u G £2(0, 1). They 
generalised the notion of the Birkhoff regular boundary conditions to this singular 
case and, for Birkhoff regular boundary conditions, found eigenvalue and eigenfunction 
asymptotics by means of the modified Priifer substitution. For the particular case of 
Dirichlet boundary conditions and the function u that is either 

(i) of bounded variation over [0, 1], or 

(ii) Lipschitz continuous on [0, 1] with exponent a G (0, 1), or 

(iii) from W?(0,1), a G [0, ±), 

the authors found several terms of asymptotic expansions for eigenvalues and eigen- 
functions. In particular, for u G W^O, 1), a G (0, |), they gave the formula 

(1.8) X n = nn — / u(t) sm(2nnt) dt + s n 

Jo 

with (s n ) G £p a for all p > 1; the remainders s n are further specified to be equal to 



;i.9) 



— / a 2 (t)cos(27mt)dt+ — / a 2 (t)dt 
-1 r t 

u(t)u(s) cos(27rnt) sin(27rns) ds dt + p r , 

Jo 



with (p n ) G t\ a . 

We also mention the papers [2,3,11-13,23,27,32], where the eigenvalue asymptotics 
were studied for several other types of singular potentials. 

In the present paper, having in mind further concrete applications to the inverse 
spectral theory, we do not aim at the utmost possible generality. Instead, we confine 
ourselves to the Dirichlet and Neumann-Dirichlet boundary conditions and potentials 
from the Sobolev space scale W^ _1 (0, 1), a G [0, 1], though the derived formula for the 
characteristic functions can be used to treat any generalised Birkhoff regular boundary 
conditions and the developed methods allow further application to potentials q = a' 
with a of bounded variation or Lipschitz continuous as in [31]. 

For the Dirichlet eigenvalues, we generalise the related results of [18] in several ways: 
a wider class of potentials is treated (no periodicity is assumed and a = is included), 
e is removed in the relation (u n ) G ^~ l ~ s for u n as in (jl.7j) (see Remark l6.2|) . and more 
terms of asymptotic expansion are given. The asymptotic formulae for derived here 
are basically the same as in [31] (see Remark 16. lj) : however, we allow the case a > | 
and simultaneously treat the Neumann-Dirichlet case. We give a special representation 
of the remainders (required for the inverse analysis) as sine Fourier coefficients of some 
functions from the W^-scale. This, however, does not yield an optimal result for the 
remainders s n of ()1.8|) in terms of the £p-scale — roughly speaking, it only implies that 
(s n ) G £l a (cf. the above mentioned results of [31]). We extract from p n of (jl.9|) one 
additional term falling into i\ a and show that the modified remainders p n form an 
^-sequence with 7 = min{3a, 1 + a}; observe that this result is incomparable to the 
inclusion (p n ) G if* proved in [31]. Also [31] gives the eigenfunction asymptotics, which 
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we do not study here (though the derived formula for the Cauchy matrix allows such 
an analysis). 

Our main result describes the asymptotics of X n and fi n in the following way. Here 
and hereafter, we denote by s n (f) and c n (f) the n-th sine and cosine Fourier coefficient 
of a function /, i.e., 

s n(f)— / f(x)sm(7rnx)dx, c n (g) = / f(x)cos(nnx)dx, n G Z + ; 



V and R will stand for the operators in L 2 (0, 1) given by 

(Vf)(x) = (1 - 2x)f(x), (Rf)(x) = f(l - x). 

Theorem 1.1. Assume that q G W / 2 a_1 (0,l) for some a G [0,1] and fix an arbitrary 
distributional primitive a G W|*(0, 1) of q. Then there exists a function a G tVf a (0, 1) 
such that 

An = 7m - s 2n {a) - s 2n (a), 

(1.10) 

/i„ = 7r(ra - |) + s 2 „-i(o r ) + S 2 n-l(<7)- 

A more detailed version of this theorem reads as follows. 

Theorem 1.2. Assume that q G W?-\0,1) for some a G [0, 1] and fix an arbitrary 
distributional primitive a G W^i®, 1) of q. Then there exists a function uj G W 2 (0, 1), 
7 := min{3o;, 1 + a}, such that 

X n — irn+ s 2n (a~) - s 2n (o-)c 2n (Va) + s 2n {u>), 

fi n = n(n - \) + s 2n _i(cr + ) - s 2n -i(a)c 2n -i(Va) + s 2n -i{uj), 



(1.11) 

where 



-l-x 



a ± (x) := ±a(x)T a\t)dt+ a(x + t)a(t) dt. 



Theorem II .21 and [31] imply the following 
Corollary 1.3. Under the assumptions (and in the notations) of Theorem \l.m\ 
(1.12) A n = 7m + s 2n (o-~) - s 2n (a)c 2n (Va) + 0(n" 7 ), n -» oo. 

In rough terms, our main result states that if the potential q belongs to the Sobolev 
space VT 2 Q_1 (0,1), then the remainders A n and Ji n are sine Fourier coefficients of a 
function from VF^O, 1). This statement agrees with the earlier known results for a = 1 
and a = mentioned above (see also [32] for the case where q is the derivative of a 
function of bounded variation, i.e., where q is a signed measure and [2,3] for Sturm- 
Liouville operators in impedance form with impedance functions from H^(0, 1)). Notice 
that Theorem 11.11 implies the related results of the paper [18], though does not imply 
the results of [31] for the operator Td and a G [0,|). However, Corollary 11.31 gives 
better than in [31] estimates for the Dirichlet eigenvalues in the uniform norm (i.e., 
in the f^-scale). Moreover, our approach is completely different from those of the 
works [18] and [31] and requires only minor effort to get the next terms in asymptotic 
eigenvalue expansions for small a. 

The organisation of the paper is as follows. In Section El we derive an equivalent 
factorised form for the differential expression l a of (jl.2j) that is more convenient for our 
purposes. We use this factorised form in the next section to derive an integral repre- 
sentation for the characteristic functions of the operators Td and Tn, and in Section 0] 
we show that the integrand in this representation possesses the desired smoothness. 
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The problem is thus reduced to finding asymptotics of zeros of certain entire functions, 
which we establish in Section and then prove the main results and some corollaries in 
Section El Several applications to the inverse spectral problem are given in Section [7| 
Finally, Appendix El provides some necessary facts from the interpolation theory. 

Throughout the paper (in addition to the above-introduced notations s n (f), c n (f) for 
sine and cosine Fourier coefficients of a function / and V and R for the multiplication 
operator by 1 — 2x and the reflection operator about x = \ respectively) W p with 
p G [l,oo) and s G R will be a shorthand notation for the Sobolev space Wp(0, 1); 
we shall also abbreviate L p (0, 1) to L p . The norm in W| is denoted by | • | s (see 
Appendix |XJ) . The space £ p consists of sequences x = (x n ) nS N with 

\x\t- ■= (^2n ps \x n \ p ^j /P < oo 

n>l 

and is a Banach space (a Hilbert space for p = 2) under the norm | • |^». For p = oo 
the norm above should be taken as \x\i^ := sup n>1 |x n |n s . 

2. Reduction of l a to the factorised form 

We recall that, for a given potential q G W% with a G [0, 1], we have defined 
the Sturm-Liouville operator Td (resp., Tn) corresponding to the differential expres- 
sion and Dirichlet (resp., Neumann-Dirichlet) boundary conditions as Tpf = 
l a (f) for / G domT D (resp., as T N f = l a (f) for / G domT N ). Here a is a fixed 
distributional primitive of q, l a is given by 

W) := -(/' - erf)' -°f' = -(-^ + °)(-^-°)f- 

and domT D and domT N are described in (jl.3j) . 

In this section, we shall derive a representation for the differential expression l a in a 
slightly different form. Roughly speaking, the claim is that the last summand (—a 2 /) 
in the above formula can be removed by changing o appropriately. More precisely, 
given r G L 2 , we denote by m T the differential expression 

•*t/)<~(sH(s-')' 

considered on the natural Z/2-domain 

domm T = {fe Wl \f-rfe Wl m T (f) G L 2 }. 

Our aim here is to show that (under some not very restrictive assumption) l a coincides 
with m T for a suitable choice of r. See also [19] for similar results on the whole axis. 

Denote by = t CT) N the quadratic form of the operator T N = T CTj n- Integration by 
parts gives that, for all / G domT N , 

t N [/] := f W)7= fW-ofr-afJ). 
Jo Jo 

It is not difficult to show (see, e.g., [5,14]) that the quadratic form L(crff' — crf'f) 
is relatively bounded with respect to the quadratic form f^df'] 2 + \f\ 2 ) with relative 
bound zero; hence Theorem VI. 1.33 of [21] implies that In is sectorial and that its 
domain is the same as for the unperturbed case a = 0, i.e., that 

domt N = {/G^ 2 1 |/(l) = 0}. 
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(As an aside we notice that, if we start from the quadratic form In and denote by 
Sn the corresponding sectorial operator, then Sn = T N ; thus the form-sum method 
and the regularisation method yield the same operator. Another consequence of this 
equality is that domT N C domt N C W%\ we also note that the inclusion domT N C W\ 
follows from the relation dom/ CT C W%, which was explained in Introduction. Finally, 
similar statements are also true for Td.) 

In the next proposition we assume that the quadratic form In is strictly accretive, i.e., 
that Re tn[f] > for all nonzero / G domtN- Since tjsi is sectorial, this situation can 
be achieved by adding to q a suitable positive constant and thus is not very restrictive. 

Proposition 2.1. Assume that a G W% , a G [0, 1], and that the quadratic form In is 
strictly accretive. Then there exists a function r G W£ such that r — a G W\ PI W 2a , 
(r — er)(0) = ; and l a = m T . Moreover, the function 

4>(x) := t(x) - a(x) + [ a 2 (t)dt, x G [0, 1], 

Jo 

belongs to with 7 = min{3a, 1 + a}. 

Proof. We shall take r in the form u'/u, where u is any function satisfying the equation 
l a (u) = and not vanishing anywhere in the interval [0, 1]. After we have proved that 
such an u exists and that r is of the required smoothness, verification of the equality 
l a = m T becomes an easy algebraic exercise (see below). 

Denote by u a solution of the equation l a (f) = satisfying the initial conditions 
w(0) = 1 and «M(0) = 0. We recall that by definition the equality l a (f) = is 
equivalent to the following first-order system: 

(2 - i} lK/ii) = (-a 2 

Since the entries of the 2x2 matrix in (j2.1j) are summable, this system enjoys the 
standard existence and uniqueness properties. In particular, the solution u with the 
stated initial conditions exists and is unique; moreover, both u and belong to Wl 
and, a posteriori, u G W\. 

We claim that u does not vanish on [0, 1]. Assume the contrary, i.e., let there exist 
x G (0, 1] such that u(x) = 0. Then integration by parts gives 

"X PX 



px fX 

= / l a (u)u = / (\u'\ 2 — auu' — au'u). 
Jo Jo 



We denote by v the function from W\ that coincides with u on [0,x] and equals zero 
on [x, 1] and observe that the above equation implies t^[v] = 0. Recall that the 
quadratic form is strictly accretive by assumption; henceforth we must have v — 0, 
which is impossible in view of the equality v(0) = u(0) = 1. The derived contradiction 
proves that u(x) 7^ for every x G [0, 1]. 

Put now <p := and r := (fi + a = u'/u. The function is in Wl and satisfies the 

equation = — (0 + cr) 2 and the initial condition 0(0) = (so that also (r — er)(0) = 0), 
and therefore 

PX PX PX 

(2.2) (j)(x) = - (f) 2 (t)dt-2 <p(t)a(t) dt - a 2 (t)dt. 

Jo Jo Jo 

It follows from Lemma fA. 31 that <p G W% and then repeated application of Lemma IA. 31 
shows that the right-hand side of (J2.2j) belongs to W 2a , so that G TV 2 2a . Next, 
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equality (|2.2jl implies that 

$(x) = - [ <p 2 (t)dt-2 I <f>(t)a(t) dt, 



and henceforth G W 7 ^ by Lemma IA.3I as claimed. 

Take now / G dom/ CT ; then f — rf = /M — <pf g W^ 1 and, using the identity 
0' = —(0 + a) 2 , we find that 

+ r) - t)/ = CfW)' - 07 - 0f + r/W - r0/ 

= (/ W)' + (0 + a) 2 / - 0/' + (0 + a)(f - a/) - (0 + a)(j>j 

= (f ll] )' + <rf' = -W). 

This shows that l a C m r . The reverse inclusion is established analogously, and, as a 
result, we get l a = m T . The proposition is proved. □ 

3. Integral representation of the characteristic functions 

Assume that a G L 2 is such that the quadratic form t^N is strictly accretive (see 
Section El for definitions). Then by Proposition 12.11 there exists r G L 2 such that 
:= r — a G Wi, 0(0) = 0, and l a = m T . Consider the differential equation m T u = X 2 u, 
i.e., 

< 3J > -(s +T )(s- T )-=^ 

It can be written first order system 

with Ui = it and u 2 = w' — ru =: ■ For any a, b G C, there exists a unique solution 
(mi,m 2 ) t of ()3.2|) subject to the initial conditions «i(0) = a, w 2 (0) = b, whence (13. 1J) 
has a unique solution w satisfying the initial conditions v(0) = a and Ut (0) = b. 

Denote by s(-, A) and c(-, A) the solutions of equation ()3.1|) obeying the initial con- 
ditions 

s(0, A) = C W (0, A) = 0, 4^(0, A) = c(0, A) = 1. 

Observe that in view of the equality 0(0) = we have wi 1] (0) = w [1] (0) - 0(O)w(O) = 
«W(0) for any u G dom/^ = domm T . Therefore the numbers ±A n coincide with the 
zeros of the entire even function s(l, •), while ±/z n coincide with those of c(l, •). In 
both cases multiplicities are taken into account (i.e., if some A 2 is an eigenvalue of Td 
of algebraic multiplicity m, then A is a zero of s(l, •) of order m, and similarly for Tn). 
The functions c(l, •) and s(l, •) are called the characteristic functions of the operators 
T N and T D respectively. 

Our aim in this section is to show that the characteristic functions c(l, •) and s(l, •) 
allow integral representations of a special form (see (|3.10|) below). We do this by 
deriving first a special integral representation for the Cauchy matrix of system (|3.2|) . 

To begin with, we notice that the matrix-valued function 

TT( \ TTt \\ ( C ^' A ) S ^' A ) 

U{x) = U(x,\):=[ w ^ 
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satisfies the initial condition £7(0) = I (with I = diag (1, 1)) and solves the equation 

(3.3) U' = (A + tJ)U, 

where 



A- A,- I ° M J- 1 ° 



In other words, {7 is the Cauchy matrix of system (|3.2jl . Since r e L 2 , equation ()3.3|1 
with the initial condition ?7 (0) = I is uniquely soluble and the solution U belongs 
to W\ entrywise. 

We shall, however, need a more explicit formula for the Cauchy matrix U. The 
standard method of variation of parameters yields the equivalent integral equation 
for U in the form 



(3.4) U(x) = e xA + e ix ~ t)A T(t)JU(t)dt, 

Jo 

where the exponent e xA can be explicitly calculated as 



X 



e 



X A _ ( cos ^ x J sm ^ x 

-AsinAa; cos Ax 



Integral equation ()3.4|) can be solved by the method of successive approximations; 
namely, with 

PX 

(3.5) U (x):=e xA and U n+1 (x) = \ e {x ~ t)A r(t)JU n (t) dt for n > 0, 

Jo 

the solution formally equals ^^L Un- Out next aim is to show that this series converges 
in a suitable topology and that the sum is indeed the Cauchy matrix. 

To this end we endow the space ^2 '■= M(2, C) of all 2 x 2 matrices with complex 
entries with the operator norm | • | of the Euclidean C 2 space and denote by W / 2 1 (^# 2 ) 
the Sobolev space of valued functions on the interval [0, 1]. 

Lemma 3.1. The series Yl^Lo Un with U n given by ()3.5j) converges in WKJI2) to the 
Cauchy matrix U. 

Proof. Bearing in mind the identity Je xA = e~ xA J and using recurrent relations f)3.5j) . 
we derive the formula 

(3.6) U n {x) = [ e^ x ~ 2 ^ A r{h) ■ ■ ■ r(t n )J n dt x . . . dt n , 

in which we have put 

]J n (x) = {t := (t u . . . ,t n ) e R n I < t n < ■ ■ ■ < t x < x}, 

n 

at) = 

Observe that < £ n (t) < x for t 6 Il„(x); thus, denoting by C the maximum of \e xA \ 
over the interval [—1,1], we get the estimate 

\U n (x)\ < C [ |r(ti)| ■ ■ ■ |r(t n )| dt l ...dt n = < ^( [ X |r| Y. 

Differentiating recurrent relations (|3.5|) . we find that 

(3.7) U' n {x) = AU n (x) + r(x)JU n ^(x) : 
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and hence, with C x := C{2\A\ 2 + 3) 1 / 2 , 

ii«/.ii^==(/(W + i^) ^(/h 2 ) (/h) • 

This estimate justifies convergence of the series Xlfclo ^ ^ n ^ ne WiG^) -topology 
to some .^-valued function V obeying the initial condition V(0) = /. Bearing in 
mind ()3.7|) and differentiating this series term-by-term, we see that V satisfies equa- 
tion ()3.3|) and thus indeed equals the Cauchy matrix U . □ 

Our next aim is to get an integral representation for U(l) of a special form. Upon 
change of variables s = £ n (t), yi = tj+i, I = 1, 2, . . . , n — 1, we recast the integral 
in ()3.6|) for x = 1 as 

U n (l)= I' e^ A r n (s)J n ds. 
Here 7~i = r(s) and, for all neff, 



o 



(3.8) r n+1 (s) = I t(s + f n (y)) r(y x ) • • • r(y n ) dy x ... dy n 
with 

(3.9) U* n (s) = {j=(y 1 ,...,y n )eR n \0<y n < y n . x < ■ ■ ■ < Vl < s + £ n (y) < !}• 

Using the Cauchy-Schwarz inequality and Fubini's theorem, we find that for every 
neN the function r n belongs to L 2 and that 

'n|o ~ 



|r n |n = / |t~„.(s) I 2 ds 



o 



1 ^ 



<7 Try I ds j \r(s + ^i(y)) r(y 1 )---r(y n - 1 )\ dy x ...dy n . 

(n-l)!y Ai;.^) 

|r(t!)| 2 • • ■ |r(t n )| 2 tftx • • • < = t \mh\ 2n 



" D! in n (i) ' " ' (n- l)!n! 



lo 



It follows that the series ^^L 1 (±l)"Tn converges in L 2 to some function r ± ; putting 
K = diag {r + ,r~}, we arrive at the desired representation for U(l): 

U(l)=e A + [ e il - 2s)A K{s)ds. 
Jo 

Spelling out the first row of this matrix equality, we get the following result. 

Theorem 3.2. Assume that o e L 2 is such that the quadratic form t CTj N is strictly 
accretive and denote by r e L 2 the function of Proposition \2. 11 for which l a = m T . 
Then the characteristic functions c(l, •) and s(l, •) of the operators and equal 

c(l, A) = cos A + / r + (s) cos A(l — 2s) ds, 
(3 - 10) • x % • xn 9^ 

smA / _ smA(l — 2s 



s(l,A) = — + J r (s) K - '-ds, 



where the L 2 -functions are defined by 

oo 

(3.11) r*= = £(±1)" 



n=l 
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with t\ = t and r n+ i given by (|H.8j) for all neN. 

In the case where a (and thus r) belongs to W 2 with some positive a the functions 
r n are also smoother. We shall establish this fact in the following section. 

4. Smoothness of the functions r n 

The derivation of the integral representations for the characteristic functions c(l, ■) 
and s(l, ■) in Section|3]only used the fact that r belongs to L 2 . If, however, the potential 
q is a distribution from W 2 ~ with some a £ (0, 1], then r £ W 7 ^ by Proposition 12.11 
and we can expect that the functions r n of (|3.8|) also have some additional smoothness. 
The aim of this section is to make this statement precise. 

Fix a natural n > 2 and consider an n-linear mapping I n : (L 2 ) n — > L 2 that acts 
according to the formula (cf. ()3.8p ) 



I n (f)(s) = / /i(s + Cn-i(y))/2(yi) • • • fn(y n -i) dyi... dy n -i, 

where f = (fi, . . ■ , f n ) € (-^2)™, the set II*_ 1 (s) is defined by (j3.9j) . and £ n _i(y) : = 
Yl?=i(~ l) l+1 yi- m particular, we see that 

r n = J„(r, . . . ,r). 
First we show that, indeed, I n maps (£2)™ into L 2 . 

Lemma 4.1. For any f = (f\, . . . , / n ) £ (L 2 ) n the function g := J n (f) belongs to L 2 
and 



l.'/l" < 7T1 IT 1/tlo- 

1 — 1 ! 



Proof. Repeating the arguments of Section El used to prove that r n £ L 2 , we find that 
(n-1)! J n; _ i(s) 



|0( S )| 2 <^-TT / \fl(s + ^-l(y))| 2 |/2(Z/l)| 2 --- l/nd/n-l)! 2 ^!...^-! 



and thus 

I <? IS < ^TTj / l/l(*l)| 2 |/2(*2)| 2 • • • Ifnitn^dt, ...dt n < f[ \fi\l 

as required. □ 

Now we give an equivalent formula for I n . Let g = I n (f) and let h be an arbitrary 
function in L 2 ; then the Z/2-scalar product (g, h)i 2 of g and h can be recast as 

(4.1) (g,h) La = [ fi(t 1 )---f n (t n )h(ut))dt 1 ...dt n . 

Jn n (i) 

Since such scalar products with h from a total set in L 2 completely determine the 
function g, we can use this formula to define the action of the mapping I n . 

The main result of this section is contained in Theorems 14.21 and 14.31 that show how 
I n acts between the spaces W 2 . 

Theorem 4.2. For a £ [0, 1], let f £ (W?) n and g := J n (f). 
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(a) For all n > 2 the function g belongs to W 2 a and there exists a positive C 
independent of f and n such that, with r : = max{0, n — 5} ; we have 

C 



l\a- 

1=1 



(b) For all n > 4 the function g belongs to and there exists a positive C 

independent of f and n such that, with r := max{0, n — 7}, we have 



C 



l9l3a -v?.W f > 



l\a- 

1=1 



Proof. Since the family {W^}, a G K, constitutes a Hilbert scale, by virtue of Theo- 
rem and Lemma f4. II it suffices to prove the theorem only for a = 1. Observe also 
that statement (a) for n = 2 holds in view of Lemma fA.4| so that we may assume that 
n > 3. 

We shall divide the proof of the case a = 1 and n > 3 into several steps and shall 
throughout denote by positive constants independent of f and n. 

Step 1. Let n > 3, f := (f u ...,/„) G {L 2 ) n and /„ G W 2 J . We shall show that in 
this case the function g belongs to W 2 l and 

„ n-1 

(4-2) Mi< — l^ l/^JJl/jlo 



for some positive constant C%. 

Let be an arbitrary test function (i.e., a C°° function with support in (0, 1)); then 
by the definition of the distributional derivative we have 

(</,0) = -(<?,<//). 

Integration by parts gives 

fn(tnW{U*)) dtn = ("I)" / fn(tnW tn {Ut)) dt n 

JO 

= (-l) ?t+1 / f n {t n )<j>(U^)) dt n 

Jo 

+ (-l) n fn(tn-M^n-2(l)) + (-l)" +1 /n(O)0 ((£n-l (t)) , 

where we have put t = (t 1; . . . , t n _i) and t = (t 1 , . . . , £ n _ 2 ). It follows now that 
-{g, 0') = (-ir +1 (/„(/!, . . . , 0) 

+ (-ir +1 /n(O)(/ n -l(/l,...,/.-l),0) 

+ (-ir(/ n _ 2 (/ 1 ,...,/ n _ 3 ,/ n - 2 ),0), 

where / n _ 2 := v4 + (/ n _ 2 , f n -i, fn) and the multilinear mapping A + is given by 
(4.3) A + (h 1 ,h 2 ,h 3 )(x) := hi(x) h 2 (y)h 3 (y) dy. 



We observe that A + acts boundedly from (L 2 ) 3 into L 2 . In fact, if h\,h 2 ,h^ G L 2 , then 
A + (h\, h 2 , h%) is a product of an L 2 -function foi and a W^-function ^h 2 h^ and thus 
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is in L 2 ; moreover, 

(4.4) |A + (/ii,/i 2 ,/i 3 )|o < |/ii|omax 

X 

Hence 



h 2 h 3 



< \hi\ \h 2 \o\h 



3 0- 



.,, (-l) n I n (fl, • • • , fn-U f' n ) + (-l) n /n(0)/„-i(/ a , . . . , f n -l) 

(4.5) n+1 

+ ( — 1)" J n _ 2 (/i, . . . , fn-3, fn-2) 



in the sense of distributions. Since the right-hand side of the above equation belongs 
to L 2 by Lemma HHJ we conclude that g G W\. 

Recall that W\ is continuously embedded into C[0, 1] and thus there is C 2 > such 
that maxj. |/| < C 2 |/|i for all / G W%- Applying Lemma I4~T1 to (|4.5jl and using the 
inequalities |/ n (0)| < C 2 \f n \i and |/„_ 2 |o < |/n~2|o|/n~i|o|/»|o (recall (jOJ)), we arrive 
at estimate (|4.2j) . 

Step 2. Let n > 3, f := (f u ...,/„) G (L 2 ) n and / x G V^ 1 . We shall show that then 
the function g belongs to H^ 1 and 



]ffli< 77=Tf l/iliIIl/: 



71 



ij\o 



i=2 



with the same G\ as in ()4.2j) . 

A direct verification shows that, with R being the reflection operator about x = |, 

^n(/l) •••)/«) — In(Rfn, ■ ■ ■ , Rfl) 

if n is even and 

-fn(/l, • • • , /n) = RIn(Rfn, • • • , -R/l) 

if n is odd. Since i? is unitary in W 7 ^ for every a G [0, 1] (for the cases a = and a = 2 
this is evident and for intermediate values follows by interpolation, see Theorem lA.lj) . 
the inclusion g G W 2 X and estimate (J4.6)) follow from the results of Step 1. 

Step 3. Let n > 3, f := (ft, ...,/„) G (L 2 ) n and / x , /„ G W^ 1 . Using relation 
and bounds (|4.2j) and (|4.6|) . we easily conclude that the function g belongs to W% and 
that 

(4-7) |p|2<-^=|/i|i|/»|xni/i|o 

1=2 

with some C3 > independent of f and r := max{0,n — 5}. The only thing to be 
justified is that for n = 3 the function h(fi) = fx — A + (fi, / 2 , fa) belongs to W\ and 

I/1I1 < C7 4 |/i|i|/ 2 |o|/ 3 |i. 
This, however, easily follows from the formula 

(A)'(z) = fxix) [ X f 3 {t)Mt) dt + /i(ar)/ a (x)/ s (x) 



showing that f[ belongs to L 2 and providing the suitable estimate of its L 2 -norm. 

Formula (|4.7|) combined with the remarks made at the beginning of the proof and 
the obvious inequality \h\ Q < \h\i completes the proof of statement (a). 
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Step 4. Let n > 4, f := (fa, ...,/„) G (L 2 ) n and fa G W^ 1 . We shall show that then 
the function g belongs to W\ and the following identity holds: 

,^ <? = —In(fu fa, fa, ■ ■ ■ , /n) — ^n-2(^.+ (/l, ^2, /s)> fa, ■ ■ ■ , fn) 

+ -^n-2(^4-(/3; /2j /l), fa, ■ ■ ■ , fn) + In-2[fl, A + (fa, fa, fa), fa, ■ ■ ■ , fa), 

where A + is the mapping of (|4.3|) and A_ is given by 

A_(h 1 ,h 2 ,h 3 )(x) := / h 2 (y)h 3 (y) dy. 



Given an arbitrary test function and integrating by parts, we get 

h / 2 W(£n(t)) tft 2 = / 2 (t 2 X(en(t)) d* 2 

^(^(^(t))^ 

+ / 2 (ti)0(e^ 2 (t)) - / 2 (t 3 )0(^2(t)), 

where t = (t 3 , t±, . . . , t n ) and t = (ti, t 4 , . . . , t n ). Substituting this relation into the 
expression for —(g,<fi') (cf. (|4.1jl ). after simple calculations we get 

- (g, 4>') = -(Ufa, fa, fa,..., fa), 0) - (I n - 2 (A+(fa, fa, fa), fa,..., fa), 0) 

+ (I n _ 2 (A_(fa, fa, fa), fa,..., fa), 0) + (J n _2 (/i, A + (fa, fa, fa), fa,..., fa), 0), 

which yields ([4.8)1 . 

Step 5. Assume that n > 4, f := (fa, . . . , fa) G (W 2 l ) n and g := I n (f). Observe that 
the multilinear transformations A + and A_ map continuously (H / 2 1 ) 3 into W 2 l , which 
can be verified by direct calculation or using Lemma IA. 31 

It follows from (j4~%|) by Step 3 of the proof that g' G W 2 2 , so that g G W 2 . The 
required norm estimate also follows from (|4.7|) and the continuity properties of A + and 
A_. Thus statement (b) of the theorem is justified for the case a = 1 and consequently, 
by the Interpolation Theorem IA.2) for all a G [0, 1]. The theorem is proved. □ 

For n = 3 we have a slightly worse result. 

Theorem 4.3. Assume that a G [0, 1], 7 := min{3a, 1 + a}, f := (fa, fa, fa) G (W?f , 
and g := /3(f). Then the function g belongs to W 2 and there is a constant C indepen- 
dent of f such that 

\g\j < Cl/iUI/^UI/sla- 

Proof. Using the definition of I3 and changing the variables via t\ = y\ — y 2 , t 2 = y 2 , 
we arrive at the representation 

9(s) = fa(s + yi- y 2 ) fa (yi) fa (y 2 ) dy 1 dy 2 

«l~s 



rs rL—s 

/ dt 2 fa(t 2 ) /i(s + ti)/ 2 (ti+t 2 )o?ti. 
Jo Jo 
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We now put 

9l (s) := J 3 (f)(s) := / dt 2 h{t 2 ) [ 2 f 2 (ti)fs(t 1 ) dh 
Jo Jo 

dt 2 f 3 (t 2 ) [ fi(t 1 )f 2 (t 1 )dt 1 
Jt 2 

and 

h(s) :=J 3 (f) := /3(f)(3) + /3(f)(3) 
and show that the function h belongs to W 2 3a and that, moreover, 

(4.9) |/i|3a<Cl|/lUI/ 2 U|/3|a 

for some C\ > independent of f. Since J 3 is a multilinear mapping, in view of 
Interpolation Theorem I A. 21 and Lemma f4. II it suffices to treat only the case a — 1. 
Assume therefore that f G (W^ 1 ) 3 . Direct calculations show that 

h'(s) = g'(s) + g[(s) = h(s) f f 2 (t)f 3 (t) dt - f / 2 (1 - + t)f 3 (t) dt 

Jo Jo 

+ / dt 2 f 3 (t 2 ) / /i(s + ti)/ 2 (ti + t 2 )dti. 
Jo jo 

Integrating by parts in the last integral, we arrive at the relation 

h\s) = -h{f 1 ,f 2 ,f 3 ), 
so that /i' G W 2 by Step 3 of the proof of Theorem 14.21 Also 

\h'\ 2 < C2I/1I1I/2I0I/2I1, 

which yields estimate ()4.9|) for a = 1. By interpolation, the results hold also for all 
intermediate a G [0, 1]. 

Consider now the function g\. Put 

J(ip u ip 2 ) := / Mt)Mt)M 
Jo 

for ipi,ip2 G L 2 ; then direct calculations show that 

9l = J(f 1,^23) ~ J(f3,Rlf>12), 

where ^23 := Ji.f2.f3) and ipu '■— J(Rf\, Rf 2 ). Applying Lemma lA.31 twice we find 
that ^12,^23 e Wi a , gi G W «+ mi ^ 2a ^ = wl and that 

min{2Q,l} min{2Q,l} — C*4 1 /l |a | /2 1 a | ^3 1 a- 

The theorem is proved. □ 

Remark 4.4. The statement of the previous theorem cannot be improved in the sense 
that the exponent 7 = min{3a, 1 + a} cannot be made larger. This follows from the 
fact that the results of Lemma IA.3I are sharp. The same statement holds also for 
r 3 = J 3 (t, r,r). 

Corollary 4.5. Assume that a G (0,1], 7 := min{3a, 1 + a}, r G W 2 and r n := 

J n (r, . . . , t). Then for every n > 3 t/ie function r n belongs to W 2 and, moreover, 

I I < C I in 
I 7 ™ 1 7 — / — r I la 

vr! 
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with some constant C > and r := max{0,n — 7}. In particular, the functions r ± in 
Theorem W . <A have the form 

r ± = ±t + r 2 + 0± 

TOt/j some W 2 -functions (f> ± . 

5. ASYMPTOTICS OF ZEROS OF SOME ENTIRE FUNCTIONS 

As we have seen in the previous sections, the eigenvalue asymptotics is completely 
determined by the asymptotics of zeros for entire functions of a special form. The main 
result of this section shows how this asymptotics can be calculated. 

Assume that / is an arbitrary function from W 2 and put 



F c (A):=cosA + / f(x) cos[A(l - 2x)) dx, 
Jo 

sin A f „, , sin[A(l — 2x)] , 
F s (\):= — + f( x )—L^ lidx. 



These are even entire functions of A; we denote by £ 2ri _i and £ 2n , n G N, zeros of F c 
and F s respectively from the set Q of ()1.4|) . We repeat every zero A 7^ according to 
its multiplicity, and if A = is a zero of F c or F s of order 2m, then we repeat it m 
times among £ 2n _i or £ 2n respectively. We shall order £ fc so that Re£ 2 n+i > Re^2n-i 
or Re^ 2n+ i = Re£ 2n _i and Im^ 2n+1 > Im^n-i and similarly for £ 2n . 

It is known (cf. [26, Ch. 1.3] and [15]) that for / G L 2 the numbers £ n have the form 

_ nH c 

for some £ 2 -sequence (^„) ne N (in particular, the remainders £ n are the Fourier coef- 
ficients of some L 2 -function). It is reasonable to expect that if / is smoother (say, 
from W2), then £ n decay faster. This is precisely what the following theorem states. 

Recall that for an L 2 -function g we have denoted by s n (g) and c n (g) its sine and cosine 
Fourier coefficients respectively, and by V : L 2 — > L 2 the operator of multiplication by 
the function (1 — 2s). 

Theorem 5.1. Assume that a G (0,1], 7 = min{3a, 1 + a}, f G Wg and that the 
numbers £ n are defined as above. Then there exists a function g G W 2 such that 

L = s n (f)-s n (f)c n (Vf) + s n (g), neN. 

In particular, (£ n ) is a sequence of sine Fourier coefficients of some function in W 2 - 

Before proceeding with the proof of the theorem, we introduce the following spaces. 
For any g G L 2 , we denote by c(g) and s(g) the sequences (c n (g))^L and (s n (g))^L 
of its cosine and sine Fourier coefficients respectively and put 

Ca := {c(g) I g G Wn, S a := {s{g) \ g G W?}, a G [0,2]. 

The lineals C a and S Q are algebraically embedded into £ 2 (Z+) and become Banach 
spaces under the induced norms 

ll c (#)llc Q := \\g\\w?, \Hg)\\s a ■= \\g\\w?- 

For any a, b G £ 2 (Z + ) we shall denote by ab the entrywise product of a and b, i.e., 
the element of £ 2 (Z + ) with entries (ab) n := a n b n . 

To establish Theorem 15 .1\ we shall essentially rely on the following three lemmata. 
The first of them ( Lemma I5.2JI is proved in Appendix El and the other two are simple 
corollaries of well known facts and thus their proofs are omitted. 
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Lemma 5.2. Suppose that a,j3<E [0, 1] and a G C a , b G Sq, ; a G Cg ; b G S^g. TTien 
aa G C a+ /3, bb G C Q+j g ; ab G S a+i g; moreover, there exists a positive constant p > 

||aa|| Ca+s < p||a|| c J|a||c a , 



(5.1) l|bb|| Ca+ , < p||b|| s J|b| 

||ab|| SQ+ ^ < p||a|| c J|b|| S; 



|9' 



Lemma 5.3. For every a G [0, 1], the operator Vf(x) = (1 — 2x)f(x) acts boundedly 
in W 2 - 

The claim follows directly from the Interpolation Theorem IA.11 
Lemma 5.4. Suppose that a G [0, 1/2). Then the Hilbert space 

oo 

H a := {a = (a n ) G £ 2 (Z+) \ a = 0, ^ n 2a \a n \ 2 < oo} 

n=l 

with norm ||a||^ a := i^2^ = in 2oi \a n \ 2 ) 1 ^ 2 coincides with the space S a , and the norms 
|| • \\ Ha and || ■ ||s Q are equivalent. 

This is a corollary of a well-known fact about Fourier transforms of spaces W 2 , see, 
e.g., [18,25,31]. 

Proof of Theorem \5.1\ Using the obvious relations 

cos£ 2 n-i = (-l) n sin£ 2n _i, 
sin£ 2n = (-l) n sin£ 2n , 
cos[6n-i(l - 2x)\ = (-l) n sin[f 2n _i(l - 2x) - (2n - l)irx)\, 
sin[6n(l - 2a;)] = (-1)" sin[f 2n (l - 2x) - 2mrx)} 

in the equalities F c (£ 2 n-i) — and F s (£ 2n ) = 0, we find that 

(5.2) sin£ n + / f(x)shx[£ n (l — 2x)—Tmx]dx = 0, neN. 

Jo 

Writing sin[£ n (l — 2x) — 7mx~\ as sin[£ n (l — 2a;)] cos(7rnx) — sin(7rnx) cos[£„(l — 2a;)], 

developing sin[£ n (l — 2x)] and cos[£ n (l — 2x)) into the Taylor series, and then changing 
summation and integration order (which is allowed in view of the absolute convergence 
of the Taylor series and the integrals), we represent ()5.2j) as 

(5.3) sin| n + 5>1)* j^yC n (V 2 ^f) - Ei-^^y^f) = °- 

k — k — 

Set £o : = 0; then, as was mentioned above, the sequence a := (£„,)nez+ belongs to 
£ 2 (Z + ), so that a G So- Define the sequence d := (d n ) n( zj 1+ through the relation 

^{2k)\ ^(2fc + l)! 
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Using Lemmata 15 .21 and 15. HI and denoting by p\ the norm of the operator V in Wg, we 
find that 



||a 2fe s(^ 2fc /)l|s Q <(pPi) 2 1|a|ir ||/||^, 
\\* 2k+l c(V 2k+1 f)h a < (ppi) 2/c+1 ||ag +1 ||/||^, 

whence d G S a . 

Equation ()5.3|) implies that sin£ n = d n for all n G Z + ; henceforth there exists n GN 
such that 

(5.4) \i n \ < 2\d n \ 

for all n > n . Fix an arbitrary number (3 G (a/3, a/2); then (3 < a/2 < 1/2 and 
d G S/3. Lemma E3 and inequalities ()5.4|) now yield the inclusion a G Sp. 
Since the sequence (sin£ n ) ne z+ can be written as 



sin a 



„2fe+l 



and since by Lemma 15.21 and the inequality 2(3 < 1 the series 

oo 

d := sina — a = ^"~^(— 1) 



„2fc+l 

\k a 



converges in S%p, we find that 

(5.5) a = d - d G S Q + S 3/3 C S a . 

By the definitions of d and d equality ()5.5|) can be recast as 

00 „2fc 00 „2fc+l 00 „2fc+l 

/c — /l — /c — 1 

Since in view of Lemmata 15.21 and 15.31 the sum 

00 „2fc 00 „2fc+l 00 „2fc+l 

falls into S 7 , we conclude that 

a-s(/) + ac(K/)GS 7l 

so that also 

[a - s(/)]c(V/) G [-ac(y/) + S,]c(Vf) C S 7 . 
This yields the desired result 

a-s(/)+s(/)c(V7)e&r, 
and the theorem is proved. □ 
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6. Proof of Theorems 11.11 fT~2l and Corollary II. Ml 

It suffices to establish only the refined asymptotics 1)1. lip . Indeed, since by Lem- 
mata and \KA\ the functions J® cr 2 (t) dt and L x a(x + t)a(t) dt belong to W 2 a and 
by Lemmata 15.21 and 15. Ml the numbers s n (a)c n (Va) are n-th sine Fourier coefficients of 
some function from W 2 a , formula ()1.1UJ1 follows from fjl.lip . 

Assume first that the potential q G W 2 ~ is such that the associated quadratic 
form In is strictly accretive. Then by Proposition 12.11 we can find a function r G W 2 
such that <p := r — a G W 2 a fl Wl and l a = m T . By Corollary 14.51 the functions 
of ()M.11)I can be represented as 

r ± = ±T + T 2 + ± 

with r 2 = / 2 (t, t ) and some H^-functions Taking into account Proposition 12. II and 
Lemma IA.4) we conclude that 

(6.1) r ± = a ± + ± , 
where 

(6.2) a ± (x) := ±a(x)T [ cr 2 (t) dt + I 2 (a,a)(x) 

Jo 

and ^ G W%. By virtue of Theorems IM.2I and 15.11 there exist functions ip^ G W% such 
that 

/n _v A n =s 2 „(r~) - s 2n (r~)c 2 „(W~) + s 2n (ip~), 

(o.M) _ 

yUn =S2n-l(^ + ) ~ S 2 „-i (t + )c 2n _l (Vt + ) + S 2 „_l(V> )• 

Since (t^ =)= cr) G VFf , equalities ()6.1|) - ()6.M|) and Lemma IB~2l yield the representation 

. t . K =S2n{<?~) ~ S 2n {(j)c 2n (V(j) + S 2n (tjj~), 

/in =S 2n -l(o- + ) - S 2 n-l(ff)c 2n -l(l / (T) + S 2 „_i(^ + ) 

with some if)^ G W^. In remains to put to := \{ip + + -R^ + + ip~ ~ to get the 

required formula. 

In a generic situation we add a suitable constant C to the potential q to get a 
potential q '.— q + C that falls into the above-considered case. (This can be done 
since the quadratic form t N is bounded below, see Section |21) Then the corresponding 
Dirichlet eigenvalues A 2 and Neumann-Dirichlet eigenvalues /i 2 have the form (J6.4)) 
with a replaced by a := a + Ct and with a^ 1 calculated as in (J6.2)) for a instead of cr. 
Since a — a G W 2 , it is easily seen that a ± — d ± E W 2 . Calculating now the integrals, 
we arrive at the relations 

A n = Tin + s 2n (cr~) - s 2n (a)c 2n (V a) + s 2n (^~), 

Ji n = 7r( n _ 1/2) + s 2n _i((T + ) - s 2n _i(cr)c 2n _i(^cr) + S 2n _i(^ + ) 

for some ^ G W 2 . Since A n = in + a n with (o n ) G £ 2 , we find that 



A„ - X n = A n - v/A2-C = + 0(A; 3 ) = ^ + - 2 
V 2A n 2vm n 2 

for some (b n ) G £ 2 , so that there exists a function x £ W^ 2 suc h that 

S2n(x) = A„ - A n , 71 G N. 
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Thus 

An = S2n(0 - S 2 n(<?)c 2 n(ycr) + S 2n ((f)~), 

with some <f>~ G W 2 . Similar arguments work for /2 n and yield the representation 

V>n = s 2n -i(o- + ) - s 2n -i(a)c 2 n-i(Va) + S 2 „_i(0 + ), 

with some <p + G W 2 . It remains to put u := \{4> + + R(j) + + </>~ — R<f>~), and the proof 
of Theorem II .21 is complete. 

Remark 6.1. Straightforward calculations show that the asymptotic formula es- 
tablished here corresponds to the asymptotic formulae (jl.8j) - (jl.9|) of the paper [31] with 
Pn = 2s 2n (o)c 2n (to) + s 2n (uj). Since (p n ) G t 2 2 a by the results of [31] and s 2n (a)c 2n (ta) 
falls into £ 2 a for a G [0, |), one concludes that (s 2n (uj)) G £ 2 a for such a. Corollary 11.31 
states that, moreover, (s 2n (u>)) G i^o f° r a ^ a e [0, 1]- 

Proof of Corollary M.cA Assume that a G W^O, 1) for some a G [0, 1] and w G W 7 ^ is 
the function of Theorem 11.21 If a G [0, 1/3], then by Lemma IA. 51 one has 

S2n(u) = 0(n~ 3a ), n — > oo, 

and, in view of (jl.llj) . 

A n = i\n + s 2n (cr~) - s 2n (a)c 2n (V a) + 0(n~ 3a ), n -> oo. 

Let now a G (1/3,1]. Since the function a G W^f satisfy then condition (vi) of 
Theorem 3.12 of [31], one has 

A n = nn- s 2n (a) + a n 
with (a n ) G l\. In view of (jl.l 1|) this yields 

S 2n (u) = -S 2n {i>) - t n , 

where 

px i-X—x 

(6.5) ^(x) := (<r~ + <r)(x) = / a 2 (t)dt+ / a(x + t)a(t) rft 



and (6„) G ^. By virtue of Lemmata IA.3I and IA.4I the function ip belongs to the 
space W 7 !". Since, moreover, ip(0) = Lemma lA.51 fsee also Remark IA.6|) yields 

the inclusion (s 2n (V'))neN £ ^2" c i- e -> ( s 2n(^))neN G ^1. Taking into account that 
7 > 1, we conclude that u>(0) = Lemma lA~5l vields now the inclusion (s 2n (uj)) G £ 2 

and the required relation ()1.12j) follows. □ 

Remark 6.2. Developing the above arguments, we can show that the function a of 
Theorem 11.11 is such that (s 2n (a)) G £ 2 a . For a G [0, |) this is shown in [31]. Observe 
that 

S2n(5 ? ) = S 2 n(^) + S 2n (a)c 2n (V a) + S 2n (V>) 

with the function uj of Theorem 11.21 and ip of ()6.5|) . For a G [|, 1] the proof of Corol- 
lary [TTSI establishes the inclusions (s 2n (u)) G £ 2 C an d ( s 2n(V0) G £2"? while 

(s 2 »c 2 „(W)) G C " % C ^ Q 

by Lemma fA. 51 

If g G W^' 1 is periodic, then we find that 

/ x co(g) c 2n (g) 

So„ (X = . 

2nV ; 2t™ 2vm 
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Inserting this into (jl.lOj) and squaring, we get 

\ 2 n = Tc 2 n 2 + co(q) - c 2n (q) + d n 

with (d n ) G ^2 Q_1 , which is to be compared to the corresponding result of [18], see (jl.7|) . 
One can also show as in [18] that this estimate is uniform in q from bounded subsets 

ofw?- 1 . 

Remark 6.3. If a = 1, then a(x) = h + f*q for some h G C, and formulae (|1.11|) give 

i r, . rl 



27m n z 27m n 



for some £ 2 -sequences (AJJ and (ju^). The choice /i = corresponds to the genuine 
Neumann boundary condition y'(0) = at x = for the operator T$, while h ^ 
produces the Robin boundary condition, y'(0) = hy(0), cf. (jl.6j) and below. 



7. Application to inverse spectral problems 

In a selfadjoint regular situation (i.e., for real- valued integrable q) the classical result 
of the inverse spectral theory states that the spectra A^ and p? n of Sturm-Liouville 
operators Tb and Tn completely determine the potential. In general the reconstruction 
algorithm is quite nontrivial and requires solvability of the so-called Gelfand-Levitan- 
Marchenko equation. In this section we shall show how Theorem 11.11 can be used in 
the inverse spectral analysis for some singular (complex-valued) potentials. 

Given the eigenvalues A^ and fi 2 n of Td and Tsr respectively, we denote by a* the 
function 

oo 

(7.1) a*(t) :=2^2(jl n sm[(2n- l)wt] - A„ sin[2vmt]) . 

n=l 

Here, as usual, X n and Jl n are defined through (jl.5j) and thus the series converges in L 2 
as soon as q G W 2 ~ l with a G [0, 1]. We also observe that, according to Theorem ll.il 
a* = a + a for the function a G W 2a of that theorem. 

Proposition 7.1. Assume that q G Wjf -1 , a G (0,1), is such that the function a* 
of (Jill) belongs to for some (3 G (0, 1], (3 > a. Then q G Wf _1 . 

Proof. We use the so-called bootstrap method. Since cr* = a + a and a G W 2a in virtue 
of Theorem ll.l[ we claim that a in fact belongs to W 2 with a 1 = min{2a,/5} > a. 
Thus the exponent a can either be taken equal to (3 or otherwise doubled. Repeating 
this procedure finitely many times, we reach the desired conclusion that a G W^, i.e., 
that q G Wi~\ ~ □ 

Developing the above arguments, we conclude that the function a* determines all 
principal singularities of the potential q. The meaning of this claim is that the functions 
a and a* share all the singularities typical for W 2 . We illustrate this issue by the 
following example. 

Assume that q is such that a(x) = J* q(t) dt has bounded variation over [0, 1] (i.e., 
da is a finite Borel measure). Slightly abusing the terminology, we shall say that the 
potential q of the corresponding Sturm-Liouville operators T ct> d and T CTj n is a finite 
Borel measure; see also [6,8,9,32] for precise definitions. 
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Proposition 7.2. Assume that q is a finite Borel measure and a* is the function 
of (J7.1)) constructed through the corresponding Dirichlet and Neumann-Dirichlet 
\j? n eigenvalues. Then the discrete parts of the measures q and da* coincide. 

i_ e 

Proof. Observe that a = J q, being of bounded variation over [0, 1], belongs to W 2 

for all e G (0, ~); henceforth a = a* — a belongs to W 2 +£ for all e G (0, ~) and 
thus is continuous. Therefore a* has the same (jump) discontinuities as a, i.e., the 
measures da* and q = da have the same discrete parts. □ 

As a final remark, we note the following. Roughly speaking, formula (|1.11|) implies 
that the Dirichlet spectrum determines the even part of the potential q, while the 
Neumann-Dirichlet spectrum the odd part of q. This explains why the unique recon- 
struction of the potential by the Dirichlet spectrum is possible if the odd part of q is 
prespecified, see [28]. 
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Appendix A. Interpolation and all that 

We recall here some facts about interpolation between W 2 spaces. For details, we 
refer the reader to [7,22,25]. 

By definition, the space W 2 coincides with L 2 and the norm | • |o in W 2 is just the 
Z/2-norm. The Sobolev space W 2 consists of all functions / in L 2 , whose distributional 
derivatives /' and /" also fall into L 2 . Being endowed with the norm 

l/|a:= {\f\l + \f\l + \f"\l) l '\ 

W 2 becomes the Hilbert space. 

Now we interpolate between W 2 2 and W 2 to get the intermediate spaces W 2 with 
norms | • | s for s G (0, 2); namely, W 2 S := [Wo 2 , W 2 ]i- a . The norms | • \ s are nondecreas- 
ing with s G [0, 2], i.e., if s < t and / G W 2 , then |/| s < |/| t . The space W 2 consists of 

L 2 functions whose distributional derivative is again in L 2 , and (|/lo + |/1o) ^ s ^ ne 
norm on W 2 l that is equivalent to | - 1 1 . Also, the Sobolev imbedding theorem implies 
that W 2 for a > 1/2 is continuously embedded into C[0, 1], the space of continuous 
functions. 

The following result [7, 22, 25] provides a powerful tool in the study of mappings 
between the spaces W 2 . 

Theorem A.l (Interpolation Theorem). Assume that < s < s lt < r < T\, and 

let T be a linear operator such that 

\Tf\ ro <C \f\ so , \Tg\ n < db| si 

for all f G W 2 ° and all g G W 2 . Put s t := (1 - t)s + tsi and r t := (1 - t)r + tr x ; 
then for every t G [0,1] the operator T acts boundedly from W 2 to with norm not 
exceeding Cl~ l C\ . 

We formulate next one result on multilinear interpolation from [7] adapted to our 
purposes. It concerns analytic scales of Banach spaces (see also [22]), and we observe 
that the scales {W 2 9+b }o & [Q ! ij are analytic for any a > and b > 0. 
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Theorem A. 2. Assume that {Eg}g e [ 0t i] and {Gg}g e [o^ are analytic scales of Banach 
spaces, neN, and that J is a multilinear mapping from (-Ei) n into G\ satisfying the 
inequalities 

n n 

\\J(s)\\g < CoH\\ 9j \\ Eo , \\J(g)\\ Gl < d J] Nk 

3=1 3=1 

for some positive constant Co, C\ and all g := (gi, g 2 , . . . , g n ) G (Ei) n . Then J can be 
uniquely extended to a multilinear mapping from (Eg) n into Go, < 9 < 1, with norm 
not exceeding Cq Cf. 

The following result is used in various places of the paper. 

Lemma A. 3. Assume that f G W 2 and g G W 2 for some a,/3 G [0,1]. Then the 
function h(x) := J x {fg) belongs to W 2 + ^ and there exists some constant C such that 

(A.l) \h\ a+f3 < C\f\ a \g\ p . 

Proof. For a = /3 = the function h is absolutely continuous and \h(x)\ < \f\o\9o\ by 
the Cauchy-Schwarz inequality, so that ()A.1|) holds with C = 1. For a = and = 1, 
we find that 

\h'\o = \fg\o < (maxb(a;)|)|/|o < C^gli, 

X 

henceforth h G W\ and (jA.lj) is satisfied for some C — C 2 > 1. To treat the case a = 1, 
(3 = 0, just interchange / and g. For a = /3 = 1 we find that h" = f'g + fg' G L 2 so 
that h G Wf with |/i| 2 < C 3 |/|i|</|i. 

Consider now a mapping M/ : L 2 — >■ L 2 given by Mfg(x) = f Q X (fg) with f E L 2 
fixed. Then, by the above, M/ acts boundedly in L 2 and W^ 1 and 

<C 2 |/| . 

By interpolation, Mf is continuous in W 2 for every /3 G [0, 1] and its norm ||Mf H^^^ 
is bounded by C 2 |/|o; in particular, 

\h\p<C 2 \f\ o \g\ . 

Analogously, for a fixed / G W 2 the operator Mf maps continuously W 2 into W 2 + 
and 

|/i|i+/3 < max{ C 2 ,C 3 } \f\i\g\p. 

The two above-displayed formulae show that M g for a fixed g G W 7 ^ is continuous as 
a mapping from L 2 into and from Vl^ 1 into W 2 1+f3 - Interpolation now yields its 
continuity as a mapping from W 2 into W 2 +f3 for an arbitrary a G [0, 1] and establishes 
inequality (jA.l|) . The lemma is proved. □ 

The next lemma establishes the required continuity of the mapping J 2 , see Section 0J 

Lemma A. 4. Assume that f G W 2 and g G W 2 for some a, f3 G [0,1]. Then the 
function 

rl—x 

h(x) := / f(x + t)g(t)dt 



belongs to W 2 +l3 and there exists some constant C such that 

\h\ a+P <C\f\ 
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Proof. The proof is completely analogous to the one of Lemma IA.3I one establishes 
first the statement for a and (3 equal to or 1, and then interpolate. The only remark 
is that for (3 = 1 one should use a representation h{x) = J^ 1 f(s)g(s — x) ds to be able 
to differentiate h. □ 

Lemma A. 5. (a) Assume that f G W 2 , a G [0, 1]; then s 2n (f) € and c 2n (f) £ £ 2 ■ 
(b) Assume that g G W% , (3 G (1,2], and g(0) = g(l); then (s 2n (g)) G i 2 . 

Proof. Part (a) follows by interpolation. Indeed, the mapping 

S : / H-> (s2n(/))neN 

is bounded from L 2 into and from W 2 into i]^. Since the spaces form an analytic 
Banach scale [22], Interpolation Theorem IA.2I yields the result about s 2n (f). Similar 
arguments justify the statement about c 2n (f). 

Part (b) requires only a slight modification. For an arbitrary s G [1, 2] we put 

Wi := {h G Wi | h(0) = h{l)}. 

Since the family {M^2 +1 }6'e[o,i] forms a Hilbert scale, by virtue of the general interpola- 
tion result from [25, Ch.l, 13.4], the family {M^2 ,+1 }eG[o,i] a ^ so ^ s a Hilbert scale. Simple 
integration by parts shows that the above operator S maps continuously W 2 into t\ 
and W| into £|. By Interpolation Theorem IA.2I the operator S maps continuously 
W 2 +1 into £ 2 +1 for all G [0, 1], and the proof is complete. □ 

Remark A. 6. Part (b) of the lemma also holds true for (3 G (|, 1) (see, e.g., [18]). 

For /^-functions / and g we denote by / * g their convolution, i.e., 

(f*9)(x) ■= f X f(x-t)g(t)dt. 



Recall also that R is the reflection operator about x = |: (Rf)(x) = /(l — x). 
Lemma A. 7. Suppose that f,g G L 2 ; then we have 

c n {f)c n {g) = c n (hx), s n (f)s n (g) = c n (h 2 ), s n {f)c n {g) = s n (h 3 ), 
where the functions hj, j = 1, 2, 3, are given by 

h x = \{R{Rf *g + f*Rg) + f*g + Rf*Rg}, 

(A.2) h 2 = l -{R(Rf *g + f*Rg)-f*g-Rf* Rg}, 

h 3 = l -{R(Rf *g-f*Rg) + f*g-Rf* Rg}. 

Proof. We shall prove only the first equality, since the other ones are treated analo- 
gously. We have 

2c n (f)c n (g) = / / f (x)g (t)[cos7m(x — t) + cos 7m(x + 1)] dxdt, 
Jo Jo 
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and simple calculations lead to 
1 /-l 

f(x)g(t) cos7m(x — t) dxdt 



o Jo 



So Oo f( s + t )9( t ) dt + J o f(t)g(s + t)dt) cosnnsds, 
f(x)g(t) cos 7m(x + £) dxdt 

= / C/ f( S ~ t } g ^ dt + f f^ 1 ~tf g ( 1 - s + t ) dt ) cosimsds. 
Taking into account the relations 

/ f(s + t)g(t)dt = R(Rf*g)(s), [ f(l-t)g(l-s + t)dt = Rf*Rg, 
Jo Jo 

we get Cn(f)c n (g) = c n (hi) with hi as stated. The lemma is proved. □ 

Proof of Lemma \5. 6 A By Lemma [A. 71 we have that c(f)c(g) = c(hi), s(/)s(p) = c(h 2 ), 
and s(f)c(g) = s(/i 3 ) for the functions hi, h 2 , and /i 3 given by (|A.2|) . 
It is easily verified that, with the operator I 2 of Section HI we have 

I 2 {f,g) = R{Rf*g). 

Recalling that R is unitary in every W 2 , we conclude by Lemma TA.4I that the functions 
hj of (jA.2j) are in W 2 +l3 as soon as / G W 2 and g G W 7 ^ 3 for some G [0, 1] and 
that, moreover, with some p > the inequality 

\hj\a+f3 < P\f\a\g\f3 

holds. Recalling the definition of the norm in S a and C a , we get the required esti- 
mates (pTff □ 



1 rl 



JO 



[1 
[2 
[3 
[4 
[5 

[6 

[7: 

[8 
[9 
[10 



References 

S. Albeverio, V. Koshmanenko, On form-sum approximations of singularly perturbed positive 
self-adjoint operators, J. Fund. Anal. 169 (1999), no. 1, 32-51. 

L.-E. Andersson, Inverse eigenvalue problems with discontinuous coefficients, Inverse Problems 
4 (1988), 353-397. 

L.-E. Andersson, Inverse eigenvalue problems for a Sturm-Liouville equation in impedance form, 
Inverse Problems 4 (1988), 929-971. 

F. V. Atkinson, W. N. Everitt, and A. Zettl, Regularization of a Sturm-Liouville problem with 
an interior singularity using quasi-derivatives, Diff. Integr. Equat. 1 (1988), no. 2, 213-221. 
Dzh.-G. Bak and A. A. Shkalikov, Multipliers in dual Sobolev spaces and Schrodinger operators 
with distribution potentials, Mat. Zametki 71 (2002), no. 5, 643-651 (in Russian); Engl, transl. 
in Math. Notes 71 (2002), no. 5-6, 587-594. 

A. Ben Amor, Ch. Remling, Direct and inverse spectral theory of one-dimensional Schrodinger 
operators with measures, Int. Equat. Oper. Theory 52(2005), 395-417. 

J. Bergh, J. Lofstrom, Interpolation spaces. An introduction, Grundlchren der Mathematischen 
Wissenschaften, Vol. 223, Springer- Verlag, Berlin-New York, 1976. 

J. F. Brasche, P. Exner, Y. A. Kuperin, and P. Seba, Schrodinger-operators with singular inter- 
actions, J. Math. Anal. Appi, 184 (1994), no. 1, 112-139. 

J. F. Brasche, R. Figari, A. Teta, Singular Schrodinger operators as limits of point interaction 
Hamiltonians, Potential Anal. 8 (1998), no. 2, 163-178. 

J. Brasche, L. Nizhnik, A generalized sum of quadratic forms, Methods Fund. Anal. Topology 8 
(2002), no. 3, 13-19. 



EIGENVALUE ASYMPTOTICS 



25 



[11] R. Carlson, R. Threadgill, C. Shubin, Sturm-Liouville eigenvalue problems with finitely many 

singularities, J. Math. Anal. Appl. 204 (1996), no. 1, 74-101. 
[12] C. F. Coleman and J. R. McLaughlin, Solution of the inverse spectral problem for an impedance 

with integrable derivative, I, Comm. Pure Appl. Math. 46 (1993), 145-184. 
[13] O. Hald, J. R. McLaughlin, Inverse problems: recovery of BV coefficients from nodes, Inverse 

Problems 14 (1998), 245-273. 
[14] R. O. Hryniv and Ya. V. Mykytyuk, ID Schrodinger operators with singular periodic potentials, 

Meth. Funct. Anal. Topol. 7 (2001), no. 4, 31-42. 
[15] R. Hryniv and Ya. Mykytyuk, Inverse spectral problems for Sturm-Liouville operators with 

singular potentials, Inverse Problems 19 (2003), 665-684. 
[16] R. O. Hryniv and Ya. V. Mykytyuk, Inverse spectral problems for Sturm-Liouville operators 

with singular potentials, II. Reconstruction by two spectra, in Functional Analysis and its Ap- 
plications, V. Kadets and W. Zelazko, eds., North-Holland Mathematics Studies, 197, 97-114, 

North-Holland Publishing Co., Amsterdam, 2004. 
[17] R. O. Hryniv and Ya. V. Mykytyuk, Inverse spectral problems for Sturm-Liouville operators with 

singular potentials, IV. Potentials in the Sobolev space scale, Proc. Edinb. Math. Soc. 49(2006) 

(to appear). 

[18] T. Kappeler and C. M6hr, Estimates for periodic and Dirichlet eigenvalues of the Schrodinger 

operator with singular potentials, J. Fund. Anal. 186 (2001), 62-91. 
[19] T. Kappeler, P. Perry, M. Shubin, and P. Topalov, The Miura map on the line, Int. Math. Res. 

Notices 2005(2005), no. 50, 3091-3133. 
[20] T. V. Karataeva, V. D. Koshmancnko, Generalized sum of operators, Mat. Zametki 66(1999), 

no. 5, 671-681 (in Russian); Engl, transl. in Math. Notes 66(1999), no. 5-6, 556-564. 
[21] T. Kato, Perturbation Theory for Linear Operators, 2nd edition, Grundlchrcn dcr Mathemati- 

schen Wissenschaften, Band 132, Springer- Verlag, Berlin-New York, 1976. 
[22] S. G. Krein, Yu. I. Petunin, and E. M. Semenov, Interpolation of linear operators, Nauka, Moskva, 

1978 (in Russian); Engl, transl.: Transl. Math. Monogr., 54 American Mathematical Society, 

Providence, R.I., 1982. 

[23] P. Kurasov and J. Larson, Spectral asymptotics for Schrodinger operators with periodic point 
interactions, J. Math. Anal. Appl. 266 (2002), no. 1, 127-148. 

[24] B. M. Levitan, Inverse Sturm-Liouville Problems, Nauka Publ., Moscow, 1984 (in Russian); Engl, 
transl: VNU Science Press, Utrecht, 1987. 

[25] J.-L. Lions and E. Magenes, Non- Homogeneous Boundary Value Problems and Applications, I, 
Springer- Verlag, Berlin-Heildelberg-New York, 1972. 

[26] V. A. Marchenko, Sturm-Liouville Operators and Their Applications, Naukova Dumka Publ., 
Kiev, 1977 (in Russian); Engl, transl.: Birkhauser Verlag, Basel, 1986. 

[27] A. McNabb, R. S. Andersscn, E. R. Lapwood, Asymptotic behavior of the eigenvalues of a Sturm- 
Liouville system with discontinuous coefficients, J. Math. Anal. Appl. 54 (1976), no. 3, 741-751. 

[28] J. Poschel and E. Trubowitz, Inverse Spectral Theory, Pure and Applied Math., Vol. 130, Aca- 
demic Press, Orlando, Florida, 1987. 

[29] A. M. Savchuk, On eigenvalues and eigenfunctions of Sturm-Liouville operators with singular 
potentials, Matem. Zametki 69 (2001), no. 2, 277-285 (in Russian); Engl, transl. in Math. Notes 
69 (2001), no. 1-2, 245-252. 

[30] A. M. Savchuk and A. A. Shkalikov, Sturm-Liouville operators with singular potentials, Matem. 
Zametki 66 (1999), no. 6, 897-912 (in Russian); Engl, transl. in Math. Notes 66 (1999), no. 5-6, 
741-753. 

[31] A. M. Savchuk and A. A. Shkalikov, Sturm-Liouville operators with distributional potentials, 
Trudy Mosk. Matem. Ob-va (Trans. Moscow Math. Soc), 64 (2003), 159-212 (in Russian). 

[32] V. V. Zhikov, On inverse Sturm-Liouville problems on a finite segment, Izv. Akad. Nauk SSSR, 
35 (1967), no. 5, 965-976 (in Russian). 

Institute for Applied Problems of Mechanics and Mathematics, 3b Naukova st., 
79601 Lviv, Ukraine and Lviv National University, 1 Universytetska st., 79602 Lviv, 
Ukraine 

E-mail address: rhryniv@iapmm.lviv.ua, yamykytyuk@yahoo.com 



